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The mean ground state occupation number and conden- 
sate fluctuations of interacting and non-interacting Bose gases 
confined in a harmonic trap are considered by using a canon- 
ical ensemble approach. To obtain the mean ground state 
occupation number and the condensate fiuctuations, an an- 
alytical description for the probability distribution function 
of the condensate is provided directly starting from the anal- 
ysis of the partition function of the system. For the ideal 
Bose gas, the probability distribution function is found to be 
a Gaussian one for the case of the harmonic trap. For the 
interacting Bose gas, using a unified approach the condensate 
fluctuations are calculated based on the lowest-order pertur- 
bation method and on Bogoliubov theory. It is found that the 
condensate fluctuations based on the lowest-order perturba- 
tion theory follow the law {S'^No) ~ A'^, while the fluctuations 
based on Bogoliubov theory behave as N'^^^ . 



I. INTRODUCTION 

The experimental achievement of Bose-Einstein con- 
densation (BEG) in dilute alkali atoms spin-polarized 
hydrogen and recently in metastable helium [p| has 
enormously stimulated the theoretical research |M-|6| on 
the ultracold bosons. Among the several intriguing ques- 
tions on the statistical properties of trapped interact- 
ing Bose gases, the problem of condensate fluctuations 
of the mean ground state occupation number 
(A^o) is of central importance. Apart from the intrin- 
sic theoretical interest, it is foreseeable that such fluc- 
tuations will become experimentally testable in the near 
future 0. On the other hand, the calculations of 
are crucial to investigate the phase collapse time of the 
condensate [|[||. 

It is well known that within a grand canonical en- 
semble the fluctuations of the condensate are given by 
= No {No + 1), implying that SNq becomes of 
order N when the temperature approaches zero. To avoid 
this sort of unphysically large condensate fluctuations, a 
canonical (or a microcanonical) ensemble has to be used 
to investigate the fluctuations of the condensate. On the 
other hand, because in the experiment the trapped atoms 
are cooled continuously from the surrounding, the system 
can be taken as being in contact with a heat bath but 



the total number of particles in the system is conserved. 
Thus it is necessary to use the canonical ensemble to in- 
vestigate the statistical properties of the trapped weakly 
interacting Bose gas. 

Within the canonical as well as the microcanonical 
ensembles, the condensate fluctuations have been stud- 
ied systematically in the case of an ideal Bose gas in 
a box jlO lJ], and in the presence of a harmonic trap 
Recently, the question of how interatomic in- 
teractions affect the condensate fluctuations has been an 
object of several theoretical investigations |p^-p7|. Idzi- 
aszek et al. [ p3[ investigated the condensate fluctuations 
of interacting Bose gases using the lowest-order pertur- 
bation theory and a two-gas model, while Giorgini et al. 

addressed this problem within a traditional particle- 
number-nonconserving Bogoliubov approach. Recently, 
Kocharovsky et al. pq | supported and extended the re- 
sults of the work of Giorgini et al. using a particle- 
number-conserving operator formalism. 

Although the condensate fluctuations are thoroughly 
investigated in Ref. p^-|26[|, to best our knowledge up 
to now an analytical description of the probability dis- 
tribution function for the interacting Bose gas directly 
from the microscopic statistics of the system has not been 
given. Note that as soon as the probability distribution 
function of the system is obtained, it is straightforward 
to get the mean ground state occupation number and the 
condensate fluctuations. The purpose of the present work 
is an attempt to provide such an analytical description 
of the probability distribution function of interacting and 
non-interacting Bose gases based on the analysis of the 
partition function of the system. 

We shall investigate in this paper the condensate fluc- 
tuations of interacting and non-interacting Bose gases 
conflned in a harmonic trap. The analytical probabil- 
ity distribution function of the condensate will be given 
directly from the partition function of the system using a 
canonical ensemble approach. For an ideal Bose gas, we 
flnd that the probability distribution of the condensate 
is a Gaussian function. In particular, our method can 
be easily extended to discuss the probability distribution 
function for a weakly interacting Bose gas. A unifled way 
is given to calculate the condensate fluctuations from the 
lowest-order perturbation theory and from Bogoliubov 
theory. We flnd that different methods of approximation 
for the interacting Bose gas give quite different predic- 
tions concerning the condensate fluctuations. We show 
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that the fluctuations based on the lowest-order perturba- 
tion theory follow the law (j^A'o) ~ N, while the fluctu- 
ations based on the Bogoliubov theory behave as iV''/^ . 

The paper is organized as follows. Sec. II is devoted to 
outline the canonical ensemble, which is developed to dis- 
cuss the probability distribution function of Bose gases. 
In Sec. Ill we investigate the condensate fluctuations of 
the ideal Bose gas confined in a harmonic trap. In Sec. 
IV the condensate fluctuations of the weakly interacting 
Bose gas are calculated based on the lowest order per- 
turbation theory. In Sec. V the condensate fluctuations 
due to collective excitations are obtained based on Bo- 
goliubov theory. Finally, Sec. VI contains a discussion 
and summary of our results. 



free energy Aq {N, Nq) and the fugacity zq of the ficti- 
tious N — Nq non-interacting bosons 



-P^Ao{N,No) = lnzo, 

Ol\o 

where the fugacity zq is determined by 

iVo = TV - V 3_ . 

^^exp[pen]zo -1 



(5) 



(6) 



We have given a simple derivation of Eqs. (||) and (g) in 
the Appendix. 

Using the free energy {N,No), the partition fimc- 
tion of the system becomes 



II. FLUCTUATIONS AND MEAN GROUND 
STATE OCCUPATION NUMBER OF THE 
CONDENSATE IN THE CANONICAL 
ENSEMBLE 



N 



Z[N]^ J2 exp[g(7V,iVo)], 



(7) 



where 



According to the canonical ensemble, the partition 
function of the system with N trapped interacting bosons 
is given by 



Z\N] = 



exp[-/3(I]„iV„e„ + (1) 



where and £„ are occupation number and energy level 
of the state n — {n^jUy^nz}, respectively. (3 — l/fc^T 
and {ux^nyjUz} are non-negative integers. Eint is the 
interaction energy of the system. For convenience, by 
separating out the ground state n = from the state 
n 7^ 0, we have 



Z[N] 



N 

E 

No=0 



q (N, No) = -p {Eo + E„,t) - pA^ {N, No) . (8) 

It is obvious that {1/Z [TV]) exp [q {N, No)] represents the 
probability finding A^o atoms in the condensate. 

To obtain the probability distribution function of the 
system, let us first investigate the largest term in the 
sum of the partition function Z [N] . Assume the number 
of the condensed atoms is A^q in the largest term of the 
partition function. The largest term exp [g (A^, A'q)] is 
determined by requiring that g^g(A^, Aq) |Aro=Arj — 0; 
i.e., 



{exp {Eo + Ernt)] Zn {N, No)}, (2) Using Eq. §) we obtain 



where Zq {N,No) stands for the partition function of a 
fictitious system comprising A^ — Ao trapped ideal non- 
condensed bosons: 



Zo (N, No) = 



E 



exp 



-f3 J2 ^n£„ 



(3) 



ln^? = /9^ {Eo+E,^t). 



(10) 



In addition, from Eq. (^) , the most probable value A"q is 
determined by 



^0 = ^ - E 

n^O exp [/?£„] (Z^) 



(11) 



Assuming Aq {N,No) is the free energy of the fictitious 
system, we have 



AoiN,No) = -fcBTlnZo(A^, A^o). 



(4) 



The calculation of the free energy ylo(A^, A^o) is non- 
trivial because there is a requirement that the number 
of non-condensed bosons is A^ — A^o in the summation 
of the partition function Zq {N,No). Using the saddle- 
point method developed by Darwin and Fowler |^^, it is 
straightforward to obtain a useful relation between the 



In the case of an ideal Bose gas, from Eq. ( p^ ) one ob- 
tains In Zq = f3eo . Thus A^q is the same as the mean 
ground state occupation number obtained by using a 
grand canonical ensemble approach. For sufficiently large 
A^, the sum X]jVo=o ^"^^ & replaced by the largest 

term, since the error omitted in doing so is statistically 
negligible. In this situation, Eq. (^l]) shows the equiv- 
alence between the canonical ensemble and the grand 
canonical ensemble for large A^. 

The other terms in the partition function (Q) will con- 
tribute to the fluctuations of the condensate, and lead to 
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the deviation of {Nq) from the most probable value A^^q. 
If iVo ^ N^, we have g^q {N, No) + 0. Assuming 



d 



dNo 



q{N,No) = a{N,No). 



(12) 



from Eqs. (||) and we obtain 



Inzo^ f3-—{Eo + E.M) + aiN,No). (13) 
dNo 

By Eqs. (|) and (|l|), we have 



No = N- 
1 



^„ exp [/3£„] exp[-(3j^ {Eo + E^^t) - a {N, No)] - 1 ' 

(14) 

Combining Eqs. (^l]) and (p^), we get the following equa- 
tion for determining a {N, No) 



No^N^o=T. 



1 



exp [Pen] exp[-/3^ {Eo + E.^)] - 1 



E 



exp [/3£„] exp[-/35^ (i;o + E^t) - a {N, No)] - 1 ' 



(15) 



Once we know Eo and Eint of the system, it is straight- 
forward to obtain a(iV, A^o) from Eq. (p^. Using 
a{N,No), one can obtain the probability distribution 
function of the system. 

From Eq. (|2|), we obtain the following result for 
q{N,No) 

I- No 

q{N,No)= a{N,No)dNo + q{N,NP). (16) 
Thus the partition function of the system becomes 

N 

Z[N]^ {exp [q {N, NP)] G {N, No)} , (17) 



No=0 



where 



G (TV, No) = exp 



No 



a {N, No) dNo 



(18) 



Assuming P {No\N) is the probability to find A^o 
atoms in the condensate, G {N, No) represents the ratio 
p(n°\^) ' relative probability to find iVo atoms 

in the condensate. From Eq. (p^, the normalized prob- 
ability distribution function is given by 



GniN,No) ^ Aexp 



No 



a {N, No) dNo 



N^ 



(19) 



where v4 is a normalization constant and is given by the 
condition A J G{N,No)dNo = 1. 

As soon as we know G {N, No), the statistical proper- 
ties of the system can be clearly described. From Eqs. 
( p7| ) and ( p^ ) one obtains the mean ground state occupa- 
tion number (iVo) and fluctuations (^d'^No) in the canon- 
ical ensemble: 



E^o=o^oexp[g(iV,iVo)] _ E^o^o ^oG (iV, iVo) 



N 



E^^^o=oexp[g(A^,iVo)] 



J:no=oG{N,No) 
(20) 



E%o=oNoGiN,No) 
Eno=oG{N,No) 



Eno=oNoG{N,No) 



J:no^oG{N,No) J 



(21) 



Starting from Eqs. (1^ and (|l|), one can calculate the 
mean ground state occupation number and fluctuations 
for ideal and interacting Bose gases. 



III. IDEAL BOSE GASES 

We now study the condensate fluctuations of the sys- 
tem with N non-interacting bosons trapped in an exter- 
nal potential. The potential is a harmonic one with the 
form 

K.*(r) = y(c.^x^+c.y+c.,V), (22) 

where m is the mass of atoms, lOx, ^y, and lUz are frequen- 
cies of the trap along three coordinate-axis directions. 
The single-particle energy level has the form 



1 



£n = { rix + - ] fiujx + \ny + - ] fiujy + I + 77 I TioJz 



(23) 



From Eq. dll]) one can get easily the most probable value 



A^Q, which reads 

np = n-n(IX 



where Tr — 



'0 huih 



is the critical temperature 
of the ideal Bose gas in the thermodynamic limit. cJ = 

1/3 

(lDx + LOy ^ LUz) /3 and ujfio — i^x^y^z) are arithmetic 



4 Ar2/3, 



(24) 
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(S'^Nq) ~ N'^/^. In this section, the Bogohubov the- 
ory will be developed based on our canonical statistics to 
discuss the condensate fluctuations originating from col- 
lective excitations. According to the Bogoliubov theory 
p9| , ^ , the total number of particles out of the conden- 
sate is given by 



•nl 



vll)fn 



(47) 



The real quantities Uni and Vni satisfy the relations 



1/2 



UnlVnl 



2f^ 



gno 



2ei 



(49) 



where /„; is the number of the collective excitations ex- 
cited in the system at the thermal equilibrium 



fn 



1 



exp [/?£ 



1 



(50) 



In addition, the energy of the collective excitations en- 
tering Eqs. ( ^8| ) and ( p9| ) is given by the dispersion law 




hujho (2n^ + 2nl + 3n + l) 



1/2 



(51) 



These phonon-like collective excitations are in excellent 
agreement with the measurement of experiments. The 
dispersion law ( pl| ) is valid if the conditions Noa/aho >> 
1 and Eni << /i are satisfied. The contribution to the 
condensate fluctuations due to these discrete low en- 
ergy modes are important because /„;, Uni+'^ni^'^niVni 
l/-\/2ri^ -I- 2nl + 3n + I at low excitation energies. 

In Eq. (^) Nni can be regarded as the effective occu- 
pation number of non-condensed atoms, while 



nl 



J nl 



(52) 



is the occupation number of the collective excitations. 
From the form of /„; one can construct the partition func- 
tion of the collective excitations in the frame of canonical 
ensemble 



Zb = ^ exp 

{nl} 



nl 



(53) 



From Eq. (|52|) Zb becomes 



{S7V„,=7V} 



7ll 



(54) 



where e'^/ — e^J (u^; -I- w^;) can be taken as an effective 
energy level of the thermal atoms. In this case Zb is 



the partition function of a fictitious boson system, which 
is composed of N non-interacting Bosons whose energy 
level is determined by s'^/ . From (^4|) the most probable 
value A^Q is given by 



exp 



-nl=0 



1 



(55) 



It is obvious that the occupation number of low n,l in 
Eq. (|5|) coincides with that of Eq. (|^. Other Nq is 
determined by 



(48) No = N - 



E — 

ni^o exp 



1 



(3 (f 



eff 



-eff 
'nl=0 



From Eqs. (|55|) and (^ we obtain 

No - 



a (TV, No) 



jnl=^0 y^-nl ~^ ^nl) fnl 



exp [-a{N,No)] - 1 
(56) 

(57) 



When getting (^7|) we have used the approximation 
fnl ~ ^bT/e^i for low energy collective excitations. Thus 
the probability distribution function of the condensate is 
given by 



Gb (iV,iVo) = ^sexp 



{No-Kf 



2 Srii/O ("ni + '"nl) fnl 



(58) 



where Ab is a normalization constant. Therefore, the 
condensate fluctuations due to the collective excitations 
reads 



Eno=oN^Gb {N,No) E'no^oNoGb {N,No) 

ENo=oGBiN,No) [ ENo=oGBiN,No) 



(59) 



Eqs. (^8|) and ( ^ ) provide the formulas for calculating 
the condensate fluctuations originating from the collec- 
tive excitations. 

Below the temperature which corresponds to the 
maximum fluctuations, we obtain the analytical result 
for the condensate fluctuations 



collective 



12C(2) 



B 



2/5 



1 „ / ma^kBTc^ 
-B ' 
2 



(60) 
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where B is a dimensionless parameter, which is the same 
as that obtained in Ref. ||2^. Note that compared with 
the result obtained by Ref. the coefhcient in ( |60| ) 
differs by a factor i. The expression ( |60| ) shows clearly 
that the condensate fluctuations due to the collective ex- 
citations are anomalous, i.e., proportional to N^^^. Note 
that Gb {N, No) is a Gaussian distribution function, the 
anomalous behavior of the condensate fluctuations comes 
from the factor 2 J2ni=£Q i^ni + fni^ which is propor- 
tional to N'^/^. 

At the critical temperature the probability distribu- 
tion is given hy Gb{T = Tc) ~ exp [— , where 



7 = 2E, 



' f^i . In this case, we obtains the 



analytical result of the condensate fluctuations 



{6^Nq) |t=t, = O.I87 = 0.18B 



ma^kBTc 



2/5 



(61) 



From ( |60[ ) and (61) we find that the behavior of the con- 
densate fluctuations based on the Bogoliubov theory is 
rather different from that of the lowest-order perturba- 
tion theory. 



VI. DISCUSSION AND CONCLUSION 

In this paper, a canonical ensemble approach has been 
developed to investigate the mean ground state occupa- 
tion number and condensate fluctuations for interacting 
and non-interacting Bose gases. Different from the con- 
ventional methods, the analytical probability distribu- 
tion function of the condensate has been obtained di- 
rectly from the partition function of the system. Based 
on the probability distribution function, we have calcu- 
lated the thermodynamic properties of the Bose gas, such 
as the condensate fraction and the fluctuations. Through 
the calculations of the probability distribution function, 
we have provided a simple method to recover the ap- 
plicability of the saddle-point method for studying the 
condensate fluctuations. In fact, the theory of the im- 
proved saddle-point method developed in this work can 
be applied straightforwardly to consider the condensate 
fluctuations in other physical systems, such as the in- 
teracting Bose gas confined in a box |4^, the interact- 
ing Bose gas in low-dimensions, etc.. The probability 
distribution function can also be used to discuss other 
interesting problems, such as the phase diffusion of the 
condensate. 

For the harmonically trapped interacting Bose gas, we 
found that different approximations for weakly interact- 
ing Bose gases give quite different theoretical predictions 
concerning the condensate fluctuations. In our opin- 
ion the lowest-order perturbation theory gives in some 
sense the condensate fluctuations due to normal thermal 
atoms, while the Bogoliubov theory gives the condensate 
fluctuations originating from the collective excitations. 



The contributions to the condensate fluctuations due to 
the collective excitations mainly come from the low en- 
ergy modes, and it is obvious that the condensate fluc- 
tuations based on the lowest-order perturbation theory 
miss the contributions coming from the collective exci- 
tations. Considering the fact that the contributions due 
to low energy thermal atoms in the lowest order pertur- 
bation theory is relatively small, the overall condensate 
fluctuations may be written as 



collective 



(62) 



where ((5^A^o)j„j and {^"^ ^0) ^^n^^f^^^ are condensate fluc- 
tuations due to the normal thermal atoms and the col- 
lective excitations, respectively. 
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APPENDIX 

In this appendix, the method of saddle-point integra- 
tion described by Darwin and Fowler [ p8| is used to in- 
vestigate the partition function of the flctitious N — A^o 
non-interacting bosons. The partition function of the fic- 
titious system is given by 



Zq (Nt) 



E 



exp 



-13 Nr^e^ 



(63) 



where Nt — N — Nq is the number of particles out of the 
condensate. 

Because of the restriction X^n/o — Nt in the sum- 
mation of Eq. (|63|), Zq (Nt) can not be explicitly eval- 
uated. To proceed we define a generating function for 
Zq (Nt) in the following manner. For any complex num- 
ber z, we take 



Go(T,z)= J2 ^""^ZoiNT) 



(64) 



Nt=0 



The generating function can be evaluated easily. The 
result of Go {T, z) is given by 



Go - Hi V 



(65) 
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To obtain Zq (Nt) we note that by definition Zq (Nt) 
is the coefficient of z^'^ in the expansion of Gq (T, z) in 
powers of z. Therefore we have 



(66) 



where the contour of integration is a closed path in the 
complex z plane about z — 0. Let g (z) be defined by 



r / M Gq (T, z) 

exp[5(z)] = ^^^+1 , 
then Zq (Nt) becomes 

Zo {Nt) ^ ^^j"^^ '^xp [g [z) 

The saddle point zq is determined by 
dg {zq) 



dzo 

From Eq. (|6^) we obtain 
d 



= 0. 



Nt 

By Eq. (pS]), one gets 



zoT^lnGo (T, Zo) - 1. 
dzo 



Nt 



E 



1 



exp [Pe^l Zq 



1 



(67) 



(68) 



(69) 



(70) 



(71) 



Noting that Eq. ([71|) is exactly the equation to determine 
the number of condensed atoms within the grand canon- 
ical ensemble, the saddle point zq can be also regarded 
as the fugacity of the fictitious N — Nq non-interacting 
bosons. 

Expanding the integrand of Eq. ( |6^ ) about z = zq, we 
have 



Zo [Nt) = ^ ^ '^^P [5 (^0) 



-(z-zo) ^g(zo) 



where 



G'^{T,zo) N^-Nt 



9zf(^°^ Go(T,zo) 
By putting z — zq — iy we obtain 



Zo (Nt) 



exp [g (zp)] 
27r 



exp 



'2dzi 



■5(^0)2/^ 



(72) 
(73) 

dy. 

(74) 



Zo (Nt) 



Go [T, Zo) 



1/2- 



(75) 



With these results the free energy Aq {N, Nq) of the fic- 
titious system is then given by 



Ao (N, No) = -ksT {In Go (T, zo) - Nt In zg 



1 



lnzo--ln [27Tg" (zo)] 



(76) 



In the case of Nt >> 1, the last two terms in Eq. ( [Tq ) 
can be omitted. Therefore 

^0 {N, No) - -keT [In Go (T, zo) - Nt In zo] . (77) 

From Eq. (|65|), we obtain the relation between 
Ao {N, No) and zo of the fictitious system: 



-l3j^Ao{N,No)^liizo. 

CiVo 



(78) 



Eqs. (|7l|) and (78) are useful relations used in the text. 



Thus we have 
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